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Active matter comprised of many self-driven units can exhibit emergent collective behaviors such 
as pattern formation and phase separation in both biological (e.g., mussel beds) and synthetic 
( e.g., colloidal swimmers) systems. While these behaviors are increasingly well understood for 
ensembles of linearly self-propelled “particles”, less is known about the collective behaviors of active 
rotating particles where energy input at the particle level gives rise to rotational particle motion. 

A recent simulation study [Nguyen et al., Phys. Rev. Lett., 2014, 112, 075701] revealed that 
active rotation can induce phase separation in mixtures of counter-rotating particles in 2D. In 
contrast to that of linearly self-propelled particles, the phase separation of counter-rotating fluids 
is accompanied by steady convective flows that originate at the fluid-fluid interface. Here, we 
investigate the influence of these flows on the coarsening dynamics of actively rotating binary liquids 
using a phenomenological, hydrodynamic model that combines a Cahn-Hilliard equation for the 
fluid composition with a Navier-Stokes equation for the fluid velocity. The effect of active rotation 
is introduced though an additional force within the Navier-Stokes equations that arises due to 
gradients in the concentrations of clockwise and counter-clockwise rotating particles. Depending on 
the strength of active rotation and that of frictional interactions with the stationary surroundings, 
we observe and explain new dynamical behaviors such as “active coarsening” via self-generated 
flows as well as the emergence of self-propelled “vortex doublets”. We confirm that many of the 
qualitative behaviors identified by the continuum model can also be found in discrete, particle-based 
simulations of actively rotating liquids. Our results highlight further opportunities for achieving 
complex dissipative structures in active materials subject to distributed actuation. 


I. INTRODUCTION 

The distributed conversion of energy into motion 
within ensembles of many self-propelled units can lead to 
complex collective behaviors operating outside the con¬ 
straints of thermodynamic equilibrium [T]. Well-studied 
examples of such active matter include migrating organ¬ 
isms |2j, the cell cytoskeleton Bill], driven granular ma¬ 
terials 00, and self-phoretic colloids [THU]. In many 
of these systems, the activity of the individual units can 
lead to phase separation and coexistence even in the ab¬ 
sence of attractive interactions. This behavior is clearly 
illustrated by simple physical models such as that of ac¬ 
tive Brownian particles (ABPs), in which hard spheres 
move at a constant speed in a direction subject to rota¬ 
tional Brownian motion [7]. ABPs are known to phase 
separate in 2D P~3HT5] and 3D [16] [17] due to a kinetic 
trapping mechanism, whereby particles incident on the 
surface of a condensed phase are trapped by other incom¬ 
ing particles m- More generally, activity-induced phase 
separation of self-propelled particles is expected when¬ 
ever the average particle velocity decreases sufficiently 
rapidly with particle density nang. This basic mech¬ 
anism is believed to underlie phase separation in such 
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disparate systems as mussel beds 03 ], bacterial colonies 

[21] . and active colloids m- 

Activity-induced phase separation has also been ob¬ 
served in systems of rotating particles, in which otherwise 
identical, gear-like disks are driven to rotate in opposite 
directions [22]. Active rotation induces effective interac¬ 
tions between the particles that can cause their segrega¬ 
tion into counter-rotating, fluid and crystalline domains 

[22] . Similar behaviors have been observed in simulations 
of spherical particles rotating within a fluid, where vis¬ 
cous shear forces couple the rotational motions of neigh¬ 
boring particles [23H25] . Experimentally, there exists sev¬ 
eral promising mechanisms by which to rotate colloidal 
components using magnetic fields [26] [27 ] . electrokinetic 
flows [28], self-phoretic motions [29H3I] , or circularly po¬ 
larized light [32]. Recent studies have shown that single 
component systems of magnetically rotated colloids |33] 
or self-rotating bacteria [53] can segregate into high and 
low density phases driven primarily by attractive dipolar 
or hydrodynamic interactions, respectively. 

Importantly, the unmixing of actively rotating par¬ 
ticles is accompanied by steady convective flows that 
originate along the interface separating the counter¬ 
rotating domains [22] . Under appropriate conditions, 
these activity-driven flows are expected to influence the 
coarsening of actively rotating fluids and may lead to 
new types of dynamically-organized structures. By con¬ 
trast, the coarsening of linearly self-propelled particles 
such as ABPs exhibit strong similarities to that of pas- 
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sive liquids [35], in which self-similar domains of size R 
grow in time as R oc t 1 / 3 . Thus, while different forms of 
microscopic activity can lead to similar mesoscale behav¬ 
iors (e.g., phase separation), these differences may con¬ 
tribute to qualitatively different dynamical behaviors at 
the macroscale (e.g., coarsening dynamics). 


discrete, particle-based simulations of actively rotating 
liquids. These results highlight opportunities for achiev¬ 
ing complex dissipative structures by directing collective 
excitations within active matter. 


II. MODEL DYNAMICS 
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FIG. 1. Phase diagram in the a(3 -plane illustrating the dif¬ 
ferent dynamical regimes (top right). Here, a measures the 
strength of active rotation, while /? measures that of frictional 
damping against the stationary surroundings (see text for de¬ 
tails) . One representative snapshot of the compositional order 
parameter p is shown for each regime. 

Here, we investigate the dynamics of liquid-liquid un¬ 
mixing in a binary fluid subject to active rotation us¬ 
ing a phenomenological, phase-field model based on the 
convective Cahn-Hilliard equation [33132 . In this de¬ 
scription, the rotational actuation of the two compo¬ 
nents in opposite directions introduces additional forces 
within the Navier-Stokes equations governing fluid mo¬ 
tion. Depending on the strength of active rotation and 
that of frictional interactions with the stationary sur¬ 
roundings, we observe three distinct dynamical regimes 
as summarized in Figure [l] For strong frictional damp¬ 
ing, coarsening of the counter-rotating domains is iden¬ 
tical to that of a passive fluid without active rotation 
(in agreement with previous particle-based simulations 
[22]). By contrast, when frictional forces are relaxed, the 
system exhibits new dynamical behaviors such as “ac¬ 
tive coarsening” driven by convective flows induced by 
the rotation of the particles as well as the emergence 
of self-propelled “vortex doublets”. We use numerical 
simulation along with scaling arguments to characterize 
the system’s dynamics within each flow regime. Addi¬ 
tionally, we show that many of the qualitative behaviors 
identified by the continuum model can also be found in 


The coarsening dynamics of binary liquids with ac¬ 
tive rotation is studied using two models: a continuum, 
hydrodynamic model and a microscopic, particle-based 
model. The former is an extension of previous phase 
field models [33133 for two phase flow that accounts for 
the active rotation of the fluid components. The micro¬ 
scopic model [22j EH] describes the Langevin dynamics 
of hard, gear-like particles, which are driven to rotate in 
opposite directions. Here, we limit our investigation to 
two-dimensional systems; however, both models can be 
readily extended to three-dimensions. 


A. Continuum Model 

In the continuum approach, we consider a binary fluid 
in which the local composition is characterized by an 
order parameter p governed by the convective Cahn- 

Hilliard equation [33133, 

| + V-M = MV 2 ft (1) 

where v is the fluid velocity, M is a mobility coefficient, 
and fi is the chemical potential. Physically, the compo¬ 
sition p(x,y,t) describes the relative amount of counter¬ 
clockwise-rotating components (p > 0) and clockwise- 
rotating components (cp < 0) at a given point in space 
and time. For simplicity, we assume the chemical poten¬ 
tial is of the form 

/i = —rep + A p 3 — K\J 2 p , (2) 

where r, A, and K are positive coefficients. These coeffi¬ 
cients determine the thickness (iF/r) 1 / 2 of the interface 
separating two phases of composition p = ±(r/A) 1 / 2 . We 
emphasize that this simple model does not attempt to 
explain the origins of phase separation driven by active 
rotation. Instead, we assume phase separation a priori 
and focus on the role of activity on the dynamics with 
which these rotating phases coarsen in time. 

To describe the activity-driven flows, we further as¬ 
sume that the fluid is incompressible, Newtonian, and 
“symmetric” such that the bulk properties of the two 
phases are equal - in particular, the density p and vis¬ 
cosity rj. Under these conditions, conservation of mass 
and momentum imply that 

V • v = 0, (3) 

p^ = -Vp + 7 ?V 2 v + pVp + V x (pr) - bv. (4) 
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In addition to the usual pressure and viscous forces 
present in the Navier-Stokes equation, equation @ in¬ 
corporates forces due to (i) capillarity [36], (ii) active 
rotation, and (iii) frictional drag, respectively. In partic¬ 
ular, we consider that the two components of the fluid 
are driven to rotate in opposite directions by a torque 
density y?r, which is proportional to the order parameter 
ip and to a constant vector r that describes the mag¬ 
nitude and direction of rotation. These local torques 
combined with spatial variations in the composition give 
rise to forces that act parallel to the interface separating 
the counter-rotating phases • In our 2D simulations, 
the fluid moves in the xy -plane with active rotation in 
the ^-direction (r = re z ). Physically, the system can 
be thought to represent an ensemble of active particles 
moving and rotating above a planar substrate as is often 
the case in experimental realizations of active matter in 
2D. To account for interactions between the particles and 
the underlying substrate, we include a frictional force in 
equation 0 characterized by a constant friction coeffi¬ 
cient b . 

At this point, it is convenient to non-dimensionalize 
the governing equations using characteristic scales for 
the interfacial thickness (LC/r) 1 / 2 , the time of unmix¬ 
ing K/Mr 2 , the equilibrium composition (r/A) 1 / 2 , and 
the chemical potential (r 3 /A) 1 / 2 . In dimensionless units, 
equations 0 and 0 reduce to 


^ + v . V ^V 2 H + /-W), (5) 

= — Vp + V 2 v + Ca _1 /iVp + aV x (ipe z ) — /3v, 

( 6 ) 


where Re = pMr /77 is a is a Reynolds number, Ca = 
MXij/K is a capillary number, and the dimensionless 
coefficients a and /3 characterize the strength of active 
rotation and frictional drag, respectively. In this paper, 
we focus exclusively on the low Reynolds number limit 
(Re —» 0) and neglect capillary forces (Ca -1 —)► 0) such 
that fluid flow is driven solely by the active rotation of 
the fluid components. 

The governing equations 0 and ([ 6 | are solved numer¬ 
ically on a square domain (L x L) with periodic bound¬ 
aries using a semi-implicit Fourier spectral method [40] 
for different values of the parameters a and /?. Initially, 
the system is prepared in a homogeneous state, in which 
the composition at each point is assigned a random value 
drawn uniformly from the interval [—0.1, 0.1]. Depending 
on the strength of active rotation a and frictional drag 
/?, this model exhibits a variety of different coarsening 
mechanisms ranging from passive, diffusive coarsening to 
active coarsening and the emergence of “vortex doublets” 
(Fig. [TJ. 


B. Microscopic Model 

To confirm the generality of active coarsening in ro¬ 
tating fluids, we study an analogous particle-based sys¬ 
tem whereby collections of hard, gear-shaped “spinners” 
are driven to rotate in opposite directions by an applied 
torque [22j [38]. Each spinner contains five circular disks 
of radius a fixed symmetrically about a central disk of 
radius 3cr. The dynamics of these composite particles 
is governed by the following Langevin equation for the 
velocity of the ith disk 

m 5 i = Fi_7Vi+F ^’ (7) 

where m is the mass of each disk, and represent 
deterministic and stochastic forces, and 7 is a frictional 
drag coefficient. The deterministic forces F^ contain both 
active and passive contributions. First, all spinners are 
driven to rotate by a constant torque = d =r with 
equal numbers rotating in each direction. Additionally, 
spinners interact both through a repulsive contact po¬ 
tential and through a short ranged attraction between 
like-rotating spinners. The latter is included to ensure 
phase separation even in the absence of active rotation 
by analogy to the continuum model, although it was not 
considered in previous works [22]. The stochastic force, 
F^ = ^ 27 fceTX(t), ensures that the system approaches 
thermal equilibrium at temperature T in the absence 
of active rotation. Langevin dynamics simulations were 
performed on graphic processing units (GPUs) with the 
HOOMD-blue [41 , 42] software package for 16, 384 spin¬ 
ners in the system. 

Although a rigorous connection between the micro¬ 
scopic and continuum models is lacking and outside the 
scope of this paper, we use order-of-magnitude reason¬ 
ing to identify dimensionless parameters a' and f3 r in the 
microscopic model that are analogous to a and /? in the 
continuum model. Specifically, a! = r/ksT measures the 
strength of active rotation relative to the thermal energy, 
whereas (3' = measures the strength of fric¬ 

tional damping. Below, all results from the microscopic 
model are presented in dimensionless form using charac¬ 
teristic scales cr, c^m/Zc^T) 1 / 2 , and ksT for length, time, 
and energy, respectively. 


III. RESULTS AND DISCUSSION 

We first use the continuum model to map out three 
qualitatively distinct parameter regimes as summarized 
in Figure [l] We discuss each regime in turn and provide 
detailed scaling arguments to explain the behaviors ob¬ 
served in the simulations. Building on insights from the 
continuum model, we reproduce many - though not all - 
of the qualitative coarsening behaviors using the micro¬ 
scopic model. 
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A. Strong Damping (f3 ^ 1) 

In the presence of strong frictional damping (/? 1), 

the coarsening dynamics of the active fluid is independent 
of the strength of active rotation (i.e., of a) and identical 
to that of a passive fluid, for which a = 0 (Fig. [2|. This 
“passive coarsening” regime has been studied extensively 
using the Cahn-Hilliard equation in the absence of fluid 
flow [43]. At short times (t <C 1), the initially homoge¬ 
neous fluid undergoes an instability characterized by a 
wavenumber k = 2 -1 / 2 , which grows in time at a rate 
1/4 until the formation of bulk domains with composi¬ 
tion ip « ±1 separated by an interfacial region of unit 
thickness. At longer times (t 1), these domains grow 
in size as R ~ t 1 / 3 due to small composition gradients 
(A (p R 4 ) which drive diffusive fluxes (j A <p/R) 
that act to grow the domains ( dR/dt ~ j) and reduce 
the curvature of the interface (Fig. [2^). Here, the do¬ 
main size R is defined as the first root of the radial pair 
correlation function [40] unless otherwise stated. 

In this regime, the active rotation of the fluid drives 
convective flows along the the interface separating the 
counter-rotating domains (Fig. [2^). The dominant terms 
of equation © are <aV x (<pe z ) ~ f3v, such that forces due 
to active rotation are everywhere balanced by frictional 
drag. As a result, all flows are directed perpendicular 
to gradients in the order parameter, and the effects of 
convective transport are negligible (i.e., v • V<p ~ 0 in 
equation (|5|). Consequently, the domain size R increases 
as R rsj independent of both a and (3 in quantitative 
agreement with passive diffusive coarsening (Fig. ©)• 

B. Weak Damping &; Weak Rotation (/3<C 1, a « 1) 

For weak frictional damping (/? 1), flows due to ac¬ 

tive rotation are no longer confined to the interface but 
rather extend into the bulk domains to influence the dy¬ 
namics of unmixing (Fig. [3|. The morphology of the 
growing domains (Fig. [3^i) is visibly different from that 
due to passive coarsening: arrays of counter-rotating vor¬ 
tices (Fig. |d) create thin filaments that break-up and 
merge with the larger domains. The characteristic do¬ 
main size R increases faster with time than expected by 
diffusive coarsening alone (Fig. [3^) . 

To understand these qualitative observations in more 
detail, consider that rotation within the bulk creates an 
interfacial stress of order a in a direction tangent to the 
interface. This active stress is balanced by the viscous 
stress U/£, where U is a characteristic velocity, and i is 
a length scale over which the velocity falls to zero. For 
small domains (R yd -1 / 2 ), velocity gradients extend 

throughout the bulk such that i ~ R and U ~ aR. As 
the domains grow larger (R /3 -1 / 2 ), the velocity de¬ 
cays exponentially with distance from the interface over 
a length yd -1 / 2 due to frictional drag; the velocity ap¬ 
proaches a constant value U a/3- 1 / 2 (Fig. S>)- ’ 

Using these estimates for the fluid velocity, we intro- 
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FIG. 2. Passive coarsening, /3 1. (a) Compositional or¬ 

der parameter ip(x, y ) at times t — 10 4 and 10 5 for parameters 
a m 10 and [3 — 10 2 ; the size of simulation cell is L — 1024. 
(b) Velocity field v(x,y) corresponding to the insets in (a). 
Arrows show the direction of fluid flow, (c) Domain size R as 
a function of time for = 10 2 and a = 0, 10, and 10 3 ; here, 
R is defined as the first zero of the radial pair correlation 
function, g(R) — 0 [4JJ- 

duce a Peclet number, Pe = £U, which characterizes 
the relative importance of convective and diffusive trans¬ 
port on the coarsening of the domains Q For small Peclet 
number (Pe < 1), activity-driven flows do not affect the 
coarsening dynamics, which is analogous to that of a pas¬ 
sive fluid. Using the above estimates for the fluid velocity, 
this condition implies that small domains, R <C a -1 / 2 , 
are unaffected by active rotation. By contrast, domains 
that grow larger than a critical size, P* ^ c^ -1 / 2 , induce 
flow velocities capable of influencing the coarsening dy¬ 
namics. This effect is illustrated in Fig. which shows 
that the domain size R(t) follows that of the passive fluid 
for R < 10 when a = 10 2 . Beyond the critical size P*, 
coarsening accelerates due to convection driven by the 
rotating fluid. In this regime, domain growth is expected 


1 In dimensional units, the Peclet number takes the more familiar 

form of Pe = tU/Mr where Mr is identified as the diffusivity. 
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FIG. 3. Active coarsening, ^ C 1 and q« 1. (a) Order 
parameter y ) at times t — 10 4 and 10 5 for parameters a = 

10 -2 and ^ = 10 —3 ; the size of simulation cell is L = 2048. 
(b) Velocity field v(x,y) corresponding to the insets in (a). 
Arrowheads show the direction of fluid flow, (c) Domain size 
R as a function of time for a = 10“ 2 and different amounts 
of frictional damping /3. 


to scale as dR/dt ~ UAcp ~ a, where Acp ~ R~ x is the 
magnitude of curvature-induced variations in composi¬ 
tion. 

Eventually, however, the rate of coarsening slows as the 
domain size R grows larger than the length £ ^ /3 -1 / 2 , 
which characterizes the decay of velocity with distance 
from the interface. Under these conditions (R /3 -1 / 2 ), 
flows are increasingly confined within a thin interfacial 
region and no longer influence the rate-limiting process 
of diffusion throughout the bulk domain. As a result, the 
domain growth returns to the diffusive scaling, R oc t 1 / 3 , 
at long times (Fig. [ 3 J 3 ). 

To summarize, “active coarsening” occurs when (i) the 
Peclet number is large, and (ii) activity-driven flows ex¬ 
tend throughout the bulk domains. These conditions are 
satisfied provided that the domain size is in the range 
a -1 ' 2 <C R <C /3 -1 / 2 . This dynamical regime is denoted 


by region II of the phase diagram in Fig. [l] 


C. Zero Damping & Weak Rotation (f3 —> 0, a <C 1) 

To better understand the “active coarsening” regime, 
we examine the limit of zero frictional damping, /? —>> 0, 
such that interfacial stresses due to active rotation are 
propagated by viscosity throughout the bulk domains 
- regardless of their size (Fig. [4|. The domain struc¬ 
ture is no longer characterized by a single length scale in 
contrast to the self-similar structures formed by “passive 
coarsening”. Instead, we observe a spectrum of different 
length scales spanning a finite range from Rmin to Rmax • 
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FIG. 4. Active coarsening, /? —> 0 and a <C 1. (a) Order 
parameter ip(x, y) at time t — 2 x 10 4 for a = 10 -2 and no 
frictional damping, /3 = 0; the size of simulation cell is L — 
2048. The right image shows a magnified view highlighting 
the two characteristic lengths, Rmin and R m ax , described in 
the text, (b) Length scales, Rmin and Rmax , as a function of 
time for a = 10 -2 and a = 10 -3 in the absence of frictional 
damping, f3 — 0. The dashed black curve corresponds to 
0.1L beyond which the finite size of the simulation domain 
becomes important; the solid black curve shows the domain 
size evolution for a passive fluid. 

The lower bound can be estimated as Rmin ~ A/C, 
where A is the total area, and C is the length of interface 
separating the bulk domains. Physically, Rmin describes 
the width of the filamentous structures that are repeat¬ 
edly drawn from the edges of the larger rotating domains. 
Because these structures are shaped by convective flows, 
their size must be sufficiently large to achieve Peclet num¬ 
bers of order unity - that is, Pe = aR ~ 1 such that 
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Rmin ~ ol 1 ' 2 . This scaling result is supported by nu¬ 
merical simulations (Fig. S>), which reveal that Rmin 
remains roughly constant throughout the coarsening pro¬ 
cess. 

The larger length scale R ma x is evaluated like R above 
as the first root of the pair correlation function. Phys¬ 
ically, the composition at two points separated by dis¬ 
tances less than Rmax are positively correlated; however, 
the strength of these correlations is considerably less than 
those observed for passive coarsening owning to hetero¬ 
geneity within these larger domains. R m ax grows roughly 
linearly in time until ca. 0.1 L, beyond which the finite 
size of the simulation cell begins to significantly influence 
the systems’ dynamics. 

In contrast to systems with frictional damping, which 
ultimately phase separate into bulk domains of arbitrary 
size, the multi-scale structures that arise in the zero- 
friction limit appear to avoid macroscopic phase sepa¬ 
ration indefinitely. Instead, active rotation continuously 
stretches and folds the growing domains in an effort to 
“mix” the fluid while it stubbornly attempts to “unmix”. 
Competition between these two processes cause the for¬ 
mation of the smaller structures of order R m in , which 
appear to persist indefinitely (barring finite size effects). 


D. Zero Damping &; Strong Rotation (f3 — > 0, a 1 ) 

As the magnitude of active rotation is increased such 
that a 1, the system transitions to a new dynami¬ 
cal regime characterized by the nucleation of localized 
vortices that move, interact, and combine within an oth¬ 
erwise homogeneous fluid (Fig. [5|. Ultimately, a single 
pair of counter-rotating vortices - a “vortex doublet” - 
emerges and propels itself autonomously throughout the 
domain, thereby mixing the fluid and preventing further 
phase separation. 

In this regime, activity-driven flows begin to shape the 
dynamics of the composition prior to the formation of 
the bulk phases. At these early times, the characteris¬ 
tic fluid velocity scales as U ~ aAcp, where Aip char¬ 
acterizes the magnitude of composition variations over a 
unit length (corresponding to the size k~ x of the fastest 
growing mode). As above, convection begins to com¬ 
pete with diffusive transport when the Peclet number is 
of order unity, Pe ^ a Aip ~ 1 . For strong rotation, 
activity-driven flows become significant even for partial 
phase separation - that is, when Aip ~ a -1 1 . 

Importantly, these convective flows have the potential 
to inhibit the further unmixing of the two fluid compo¬ 
nents. In the absence of active rotation, fluid unmixing 
proceeds exponentially as Aip oc exp(£/4). To inhibit 
phase separation, the shear rate in the fluid must ex¬ 
ceed the rate of unmixing. Partial phase separation re¬ 
sults in activity-driven flows with shear rates of order 
a Aip. Thus, when the extent of unmixing reaches a crit¬ 
ical value - namely, Aip > <a _1 - the resulting flows will 
act to oppose further unmixing. 
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FIG. 5. Vortex doublet, /5 — > 0 and a ^$> 1. (a) Order 
parameter ip(x,y) at times t — 20, 48, 150, and 300 for pa¬ 
rameters a = 50 and /3 = 0; the size of simulation cell is 
L = 256. Number of vortices as a function of time for ft = 0 
and a — 25 and 50. Here, a vortex is identified as a localized 
region in which the composition is \<p\ > 0.3. 


Nevertheless, the spatial heterogeneity of the shear 
field allows for the nucleation of small vortices in lo¬ 
cations with lower shear rates. As vortices form, they 
create regions of low shear (but high vorticity) in their 
interior that allow for further phase separation. At the 
same time, these vortices induce high shear rates in the 
surrounding fluid, which inhibits unmixing therein. As a 
vortex strengthens, the composition in its interior ap¬ 
proaches ip ~ ±1, while that of the exterior remains 
largely homogeneous with \ip\ 1. 

Following the initial nucleation phase, vortices move 
and deform in the swirling flows induced by their neigh¬ 
bors. Some are destroyed by strong shear due to larger 
neighbors; others grow and merge to form larger and 
more powerful vortices. In this way, the number of vor- 
ticies decreases in time until only two large, counter¬ 
rotating vortices remain (Fig. [5]). Together, these vor¬ 
tices form a stable, self-propelled vortex doublet of size 
R that swims about the domain with velocity U ~ aR. 
The doublet creates a velocity disturbance that decays as 
aR 3 /r 2 with distance r from its center. Consequently, 
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a single vortex doublet can create shear rates of order 
unity (that necessary for mixing) at distances of Ra 1 / 3 . 
This result is consistent with the simulation results shown 
in Figure [5j in which a single vortex doublet effectively 
mixes a region ca. four times as large as itself with 
a = 50. 


(a) 


E. Results of the microscopic model 



The key insight suggested by the continuum model 
is that phase separation in actively rotating liquids can 
drive convective flows that feedback into the system and 
direct the dynamic evolution of the growing phases. This 
result is further supported by microscopic simulations of 
actively rotating particles (Fig. |6|. In particular, we 
studied the collective dynamics of 16, 384 spinners sub¬ 
ject to moderate driving torques (a! = 0.25) and differ¬ 
ent levels of frictional damping. Under these conditions, 
spinners unmix to form domains of like-rotating particles 
(Fig. [6 ^l, b) that grow steadily in time. Here, the size R' 
of the growing domains is quantified as the first root of 
the integral of the pair correlation function. Note that 
this measure is different from that used in the contin¬ 
uum model and is chosen for its decreased sensitivity to 
statistical fluctuations. 

For relatively strong damping (/T = 1), the domain 
size increases as R' oc t' 1 / 3 (Fig. pa,c) in agreement with 
the continuum model (Fig. [2| and with previous micro¬ 
scopic results [22]. By contrast, as the damping param¬ 
eter is decreased, convective flows emerge and accelerate 
the rate of coarsening (Fig. Consistent with the 

continuum model, the growth exponent increases from 
1/3 toward 1 as the damping parameter /3' is reduced. 
The exploration of smaller damping parameters and/or 
longer simulations times in the microscopic model was 
found to be computationally prohibitive. Therefore, it is 
unclear if the microscopic model will approach a scaling 
exponent of ca. 1 in the limit as ft —>• 0 or if it will return 
to a scaling exponent of 1/3 in the limit of long times. 

To explore the possibility of “vortex doublets” in the 
microscopic model, we increased the driving torque to 
a' = 2.5 under conditions of weak damping (ft = 0.01). 
Consistent with the continuum model, the increased rota¬ 
tional activity of the particles was sufficient to inhibit the 
unmixing of the spinners; however, we did not observe the 
nucleation of localized vortices from homogeneous initial 
conditions. Furthermore, we applied the same driving 
torque to an initially phase-separated system under con¬ 
ditions of weak damping {a! = 10 and ft = 0.01). The 
active rotation of the particles resulted in their complete 
mixing, which suggests that the absence of vortex dou¬ 
blets in the microscopic model is not the result of a nu¬ 
cleation barrier. 

The discrepancies between the two models under con¬ 
ditions of high torque likely arise from a failure of the 
continuum model to account for the microscopic effects 
of active rotation on fluid-fluid phase separation. In the 



FIG. 6. Microscopic model. (a,b) Representative snap¬ 
shots of a 50:50 mixture of 16,384 spinners driven to rotate 
in clockwise (orange) and counterclockwise (purple) directions 
at time t'm 10 5 . The visualization of the particles by their 
voronoi tessellation is shown in the right images along with 
the fluid flows. The image in (a) show “passive coarsening” 
for parameters a'm 0.25 and ft — 1 corresponding to strong 
frictional damping. Image in (b) show “active coarsening” for 
parameters a — 0.25 and ft = 0.01 corresponding to weak 
damping. White arrows show the direction of fluid flows, (c) 
Domain size R' as a function of time for ol = 0.25 and differ¬ 
ent amounts of frictional damping ft — 1, 0.05, and 0.01; this 
plot is analogous to that in Fig. [3p. Here, R' is defined as the 
first zero of the integral of the pair correlation function. 


microscopic model, the driving torque a' sets an energy 
scale, which must be significantly less than that of at¬ 
tractive interparticle interactions to achieve microscopic 
phase separation. In the high-torque simulations, these 
two energy scales are comparable such that active ro¬ 
tation prohibits any and all phase separation. In such 
systems, the strong torques required for the formation 
of vortex structures lead instead to the complete mix¬ 
ing of the binary fluid. Additional theoretical work is 
required to account for the microscopic details of the 


















particle-based simulations (e.g., particle shape) within 
an accurate hydrodynamic description. Nevertheless, it 
remains likely that the vortex structures observed in the 
hydrodynamic model could also be realized in analogous 
microscopic models provided the driving force for phase 
separation is sufficiently strong. 

IV. CONCLUSIONS 

To summarize, we presented a phenomenological, con¬ 
tinuum model for studying the dynamics of phase sep¬ 
aration in binary liquids with active rotation. Convec¬ 
tive flows induced by the fluid activity result in acceler¬ 
ated coarsening as compared to spinodal decomposition 
in passive liquids. The transition from passive to active 
coarsening is determined primarily by the strength of ac¬ 
tive rotation and that of frictional interactions with the 
stationary surroundings. In addition to active coarsen¬ 
ing, the continuum model also predicts the formation of 
self-propelled vortex doublets under conditions of strong 
rotation and weak frictional damping. These dissipative 
structures emerge spontaneously from the competition 


between fluid mixing via active rotation and fluid unmix¬ 
ing due to interparticle interactions. Many of the trends 
observed in the continuum model such as accelerated 
coarsening are also reproduced by microscopic kinetic 
simulations of counter-rotating particle mixtures. Fur¬ 
ther work is needed to develop a more rigorous connection 
between such microscopic models and the continuum hy¬ 
drodynamics of actively rotating fluids. We are currently 
developing experimental models of counter-rotating par¬ 
ticle mixtures to explore and elucidate the collective dy¬ 
namics that emerge from active rotation. Such nonequi- 
librium materials have the potential to achieve life-like 
properties such the abilities to adapt, reconfigure, and 
repair dynamically in response to environmental stimuli. 
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